Comment on 

"Gravitationally Induced Neutrino-Oscillation Phases" 



Tanmoy Bhattacharya 1 , Salman Habib 2 , and Emil Mottola 1 



1 Theoretical Division 

T-8, MS B285 
Los Alamos National Laboratory 
Los Alamos, NM 87545 

2 Theoretical Division 
T-DOT, MS B288 

Los Alamos National Laboratory 
Los Alamos, NM 87545 



Abstract 

We critically examine the recent claim of a "new effect" of gravitationally 
induced quantum mechanical phases in neutrino oscillations. A straightfor- 
ward exercise in the Schwarzschild coordinates appropriate to a spherically 
symmetric non-rotating star shows that, although there is a general relativis- 
tic effect of the star's gravity on neutrino oscillations, it is not of the form 
claimed, and is too small to be measured. 



In a recent communication M, Ahluwalia and Burgard claim to have 
discovered a "new effect from an hitherto unexplored interplay of gravitation 
and the principle of the linear superposition of quantum mechanics." In fact, 
the calculation of the quantum mechanical phase of a particle propagating 
in the geometry of a collapsed star appears in several textbooks on general 
relativity 0. More importantly, the claimed results in equations (6) to (8) 
of Ref. [ l] appear to be at variance with the standard treatment, which we 
review in this Comment. The gravitational correction to the phase turns 
out to be of order 10~ 9 for neutrinos in the eV mass range and is therefore 
completely negligible for situations of astrophysical interest. 

In the geometrical optics limit the quantum mechanical phase accumu- 
lated by a particle propagating from point A to point B in the gravitational 
field described by the metric g^ u is given by the classical action of the particle, 
divided by h, namely, 
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$ = -/ mds = - Pudx 11 = - (-Edt + pidx 1 ) (1) 
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where p^ is the four momentum conjugate to x^\ 

Pn = mg^^-, (2) 

and ds is an element of proper length of the particle's wordline. The integrand 
of Eqn. ([I]) is obviously an invariant under coordinate transformations. How- 
ever, the form of the integral will depend on the labeling of the end-points 
A and B. Equation ([!]) is the same as the Eqn. (4) of Ref. |lj with which 
those authors begin. 

The authors of Ref. address the radial propagation of relativistic 
neutrinos in the potential of a spherically symmetric non-rotating star which 
is described by the Schwarzschild line element 

ds * = -(l- ^) dt 2 + ^) ' dr 2 + r 2 (d6 2 + sin 2 9 d</> 2 ) . (3) 

We note that the semiclassical phase for radial motion in a spherically sym- 
metric background does not depend on the spin of the particle, as can be 
verified by explicit calculation using the spin connection in the Dirac equa- 
tion in this background [@ . Hence Eqn. ([!]) applies equally well to neutrinos 
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as to scalar particles in the case of radial motion. Because the spacetime 
has a time-like Killing vector, d/dt, the momentum conjugate to it is time 
independent: 

/ 2GM\ dt 
E = —p t — mil I — = constant . (4) 

The value of this constant E is the asymptotic energy of the neutrino at 
r = oo. For radial motion, the mass shell constraint reads 

7 / 2GM\~ 1 „ / 2GM\ , , 
9' P»Pu + m 2 = = - ( 1 — J E + f 1 —J p 2 r + m 2 , (5) 

from which we obtain 



r 

Making use of the definitions 



, 2GM\ / 9 _ 2GMm 2 
p r (l )=\E 2 -m 2 + . (6) 



Pr = m[l-—) -, (7) 



and (HI), we can write 



dr \ r J p r 

We regard the Schwarzschild radial coordinates and as fixed and express 
the phase $ in those coordinates as follows: 
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This is just the standard result 0. In the weak field approximation, we 
expand this result to first order in G, to obtain 



m 2 . , GMm 4 , fr B \ , . 

$ ~ (r B - r A ) H ^ In — + • • • . (10) 

fiV£ 2 -m 2 ft(£ 2 -m 2 ) 3/2 W 1 ; 

The precise application of the above formula depends on the physical situa- 
tion at hand, specifically, on what variables are to be held fixed in a particular 
interference experiment. 

Let us consider the case of neutrinos produced at fixed asymptotic energy 
E in a weak flavor eigenstate that is a linear superposition of mass eigen- 
states, mi and mi- Since the energy is fixed but the masses are different, if 
interference is to be observed at the same final spacetime point (r^tg), the 
relevant components of the wave function could not both have started from 
the same initial spacetime point (ta, ^a) m the semiclassical approximation. 
Instead the lighter mass (hence faster moving) component must either have 
started at the same time from a spatial location r < ta, or (what is equiva- 
lent) started from the same location ta at a later time t A + At. Hence, there 
is already an initial phase difference between the two mass components due 
to this time lag, even before the transport from r A to r B which leads to the 
phase $ in (0). The additional initial phase difference may be taken into 
account most conveniently from the second point of view, i.e. by treating 
the spatial coordinates r A and r B as fixed and the time of transit, 

dt= -dr (11) 



A Jr A dr 

as the dependent variable through Eqns. @ and (|j). The difference of this 
time of transit between the two mass eigenstates, multiplied by E is precisely 
the additional phase, EAt which we must add to A$ to obtain the correct 
relative phase between the two mass components of the wave function which 
interfere at (r B , t B ) with fixed energy E. We note that the COW experiment 
0] may be treated by similar reasoning and that many other scenarios for 
neutrino oscillations which may be envisaged lead to the same result. 
Hence we are led to compute instead of the quantity, 

rB rB 

$ r = $ + £ / dt = p r dr (12) 
J A J a 
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This $ r may be calculated just as easily as the full $: 



1 f s JE 2 -m 2 + 
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In the weak field expansion this becomes 
(r B -r A ) + 2GMln 
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h\ / E 2 -m 2 \r A , 



(14) 

The relative phase difference A$ r between two mass eigenstates of relativistic 
neutrinos (E 2 ^> m 2 ) with mass squared difference of Am 2 created at point 
ta and interfering at point r# is then 

(15) 

where c has been restored to facilitate numerical calculations. We note that 
in the relativistic limit this result is precisely minus half of the equivalent 
quantity computed from the full phase $. The first term (A$°) in Eqn. (|15|) 
is the standard flat space result, well known in both neutrino and strangeness 
oscillations. The leading order correction to this familiar result has cancelled 
in the relativistic limit and we are left only with the latter higher order terms 
in Eqn. ([15]). The second term is the special relativistic correction to the 
phase which is usually neglected for light neutrinos, and the last term is 
the effect of the gravitational field of the star in static Schwarzschild coordi- 
nates which enters only at the same higher order in 1/E 3 . Numerically its 
magnitude is equal to 

which is completely negligible in typical astrophysical applications. 

We note more generally that for the case of radial motion in coordinates 
such that g rt = 0, the expression for the phase, Eqn. ([[]), can be written as 



m 2 r B dry, 
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where 




V-gttdt 
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and 




(18) 



which is a statement of the equivalence principle. Locally there are no ob- 
servable effects of the gravitational potential. The redshifted energy -Ei OC ai 
is not a constant of motion and cannot be pulled outside of integrals, which 
accounts for the appearance of (l/E locaX ) av = (1/T AB ) £f dt\ OC3 \/ E locai in the 
above expression. 

The authors of Ref. [I] claim to find an effect on the phase, first order in 



However no derivation is given to support this claim and their quantity E 
is never defined. If E is the constant of the motion defined in Eqn. (01), 
then this claim disagrees with the standard result (|T5| ) rederived here, and 
is therefore incorrect. On the other hand, if E is to be identified with our 
£i oca i, then it cannot be pulled out of the integral, and (|T9|) is incorrect for 
that reason. 

1 Actually the phase quoted in Ref. Q Eqn. (6) differs from the one obtained by us by 
a factor of two, even in fiat space, because A$ r enters interference probabilities through 
cos (A$ r ) whereas their Eqn. (6) contains sin 2 of a phase angle which should be A<I> r /2 
by the half-angle formula. This discrepancy may have arisen from a confusion between $ 
and <£v which also differ by a factor of — ^ , as noted previously. 



GQ, equal to (Eqn. (8) of Ref. 0) 




(19) 
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There is a sense in which the first (A$°) term of Eqn. ( |T5D has a contri- 
bution similar in form to (1l9[). If we define 



- 1 f A Edr 

-C/local = / i , (20) 

tb - r A Jb ,/l _ 2GM 



then in the weak field limit 



EGM , /r B \ 
+ In — 21 



and we can write 



A ^ n Am 2 c 3 , v Am 2 GMc . fr B \ 

Zn&iacgX ^local \^A/ 

which is of the form reported in Ref. M. However, to interpret this rewriting 
of the standard result in different variables as a "new" gravitational effect 
would be very misleading. The quantity -Bi OC ai refers to the energy measured 
by local observers at fixed r and differs from the asymptotic E precisely 
because of the well known gravitational redshift effect of general relativity. 
Since all measuring rods and clocks are subject to the same redshift, there 
are no physical consequences of this local redshift effect on the observable 
physics of neutrino oscillations (for example, if along with the local energy 
one uses the proper, rather than the coordinate, length, then the "effect" in 
( |22| ) disappears). In particular, any effect (s) of neutrino oscillations on energy 
transport and heating in supernova explosions are quite indifferent to such 
local redefinitions of length, time and energy scales. Provided all calculations 
are done in a relativistically covariant framework, local redshift effects are 
accounted for automatically and there are no observable consequences for 
supernova evolution to be deduced from the decomposition in (|22|). Of course, 
if one does not use a relativistically covariant framework in the calculations, 
the error made will be precisely of the order of the second term in (p2|). 

The only indication of the basis for the claim in Ref. [|TJ is a reference 
to a paper by Stodolsky [[|. However, as Stodolsky himself notes, the split 
between "flat" and "curved" space effects in equation (2.3) of his paper is 
coordinate dependent. Hence there is no invariant meaning to the splitting of 
the phase into these two pieces, and such a splitting is completely misleading 
for the present application, just as is the splitting of A$° into the two pieces 
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in (|22| ) above. In addition, the time component of the quantity Stodolsky 
calls the "usual four-momentum of special relativity" is not a constant of 
motion in the present application and cannot be removed from integrals over 
r. Since Stodolsky starts with precisely the same phase $ of Eqn. ([I]) the 
sum of his two pieces is precisely equal to the same result @ rederived here, 
as may be checked directly from the definitions in Ref. ||. 

Finally we note that the gravitational effect which we have computed here 
in ([16]) has a different dependence on the neutrino masses and energy from 
the flat space result, A$°, and hence it cannot be absorbed into A$° by a 
coordinate transformation. 
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